rm(list=ls()) y0 <- 10 T <- 100 time <-1:T g <- 0.05 y <- rep(0,T) y_L <- rep(0,T) y_g <- rep(0,T)